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A TWO-SPECIES COMPETITION MODEL ON Z D 


GEORGE KORDZAKHIA AND STEVEN P. LALLEY 


Abstract. We consider a two-type stochastic competition model on 
the integer lattice Z d . The model describes the space evolution of two 
“species” competing for territory along their boundaries. Each site of 
the space may contain only one representative (also referred to as a 
particle) of either type. The spread mechanism for both species is the 
same: each particle produces offspring independently of other particles 
and can place them only at the neighboring sites that are either un¬ 
occupied, or occupied by particles of the opposite type. In the second 
case, the old particle is killed by the newborn. The rate of birth for 
each particle is equal to the number of neighboring sites available for 
expansion. The main problem we address concerns the possibility of 
the long-term coexistence of the two species. We have shown that if we 
start the process with finitely many representatives of each type, then, 
under the assumption that the limit set in the corresponding first pas¬ 
sage percolation model is uniformly curved, there is positive probability 
of coexistence. 

Key words: coexistence, first passage percolation, shape theorem. 

1. Introduction 

The voter model is an interacting particle system in which individuals 
(particles) of two species, Red and Blue, compete for “territory” on a (locally 
finite) graph. At each time t > 0, every vertex (site) of the graph is occupied 
by a single particle, either Red or Blue. At any time, a particle of color i at 
a vertex x may spontaneously die, at rate equal to the degree of x, and be 
replaced by a clone of a randomly chosen neighbor. Thus, a vertex of color 
i spontaneously flips to the opposite color j at rate equal to the number 
of neighboring vertices of color j. See 7] for a formal construction of this 
process and an exposition of its basic properties. 

The Richardson model [S; was introduced as a model for the spatial spread 
of a population in a favorable environment. The environment is once again 
a locally finite graph. At any time a vertex may be occupied by at most 
one particle (some vertices may be unoccupied); all particles are of the same 
species. Once occupied, a vertex remains occupied forever. Each unoccupied 


Date : February 2, 2008. 


1 



2 


GEORGE KORDZAKHIA AND STEVEN P. LALLEY 


vertex is spontaneously occupied at instantaneous rate equal to the number 
of occupied neighbors. 

In this paper we study a hybrid of the voter and Richardson models on 
the integer lattice Z rf , which we dub the two-species competition model , or 
simply the competition model. The dynamics are as in the voter model, but 
unlike the voter model, vertices may be unoccupied. An unoccupied vertex 
is colonized at rate equal to the number of occupied neighbors, as in the 
Richardson model; at the instant of first colonization, the vertex flips to 
the color of a randomly chosen occupied neighbor. Once occupied, a vertex 
remains occupied forever, but its color may flip, as in the voter model: the 
flip rate is equal to the number of neighbors occupied by particles of the 
opposite color. The state of the system at any time t is given by the pair 
( R(t),B(t )), where R(t ) and B(t) denote the set of sites occupied by Red 
and Blue particles, respectively. Note that the set R{t) U B(t) of occupied 
sites evolves precisely as in the Richardson model, and so the growth of this 
set is governed by the same Shape Theorem (see section 11.11 below) as is the 
Richardson model. 

Our primary interest is in the possibility of long-term coexistence of the 
two species, given initial conditions in which only finitely many vertices 
are occupied (with at least one vertex of each color). It is clear that at 
least one of the two species must survive, and that for any nondegenerate 
finite initial configuration of colors there is positive probability that Red 
survives and positive probability that Blue survives. However, it is not at 
all obvious (except perhaps in the case where the ambient graph on which 
the competition takes place is the integer lattice Z - see section |3 below) 
that the event of mutual survival has positive probability. Our main result 
concerns the competition model on the graph Z^. Say that a compact, 
convex set S with boundary dS is uniformly curved if there exists g < oo 
such that for every point z G dS there is a ball of radius g with z on its 
surface that contains S. 

Theorem 1. If the limit shape S for the Richardson model is uniformly 
curved, then for any nondegenerate initial finite configuration the event of 
mutual survival of the two species has positive probability. 

The proof will be carried out in sections El El below. Theorem 0 is by no 
means a complete solution to the coexistence problem, because it remains 
unknown whether the limit shape S for the Richardson model is uniformly 
curved, or even if its boundary dS is strictly convex. Nevertheless, simula¬ 
tions give every indication that it is, and jBJ suggests a possible explanation 
of what lies behind the strict convexity of dS. 

The two-species complete model is superficially similar to the two-type 
Richardson model studied by Haggstrom and Pemantle [3j, but differs in that 
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it allows displacement of colors on occupied sites: in the two-type Richardson 
model, once a vertex is occupied by a particle (either red or blue) it remains 
occupied by that color forever. The main result of |lj is similar to Theoremd 
but requires no hypothesis about the Richardson shape: it states that mutual 
unbounded growth has positive probability. Because no displacements are 
allowed, the behavior of the two-type Richardson model is very closely tied 
up with the first-passage percolation process with exponential edge passage 
times. The two-species competition model is also closely related to first- 
passage percolation, but the connection is less direct, because the possibility 
of displacements implies that not only the first passages across edges play a 
role in the evolution. 

We have run simulations of the competition model with initial configu¬ 
ration R( 0) = {(0,0)} and £?(0) = {(1,0)}. Figured shows two snap shots 
of the same realization of the process taken at the times when the region 
occupied by both types, R(t) U B(t), hits the boundary of the rectangles 
[-300, 300] x [-300, 300] and [-400, 400] x [-400, 400] respectively. Ob¬ 
serve that the overall shape of the red and the blue clusters did not change 
significantly. We believe that the shape of the regions occupied by the red 
and blue types stabilizes as times goes to infinity. 

For any subset Z C define 

Z = {x G : dist(x, Z) < 1/2}, 

where dist denotes distance in the L 00 —norm on M d . For any subset Z C 
and any scalar s > 0, let Z/s = {y/s : y E Z}. 

Conjecture 1. There exist random sets R and B such that with probability 


one 


(1) 

lim R(t)/t = R, 

t —»oo 

(2) 

lim B(t)/t = B, 

t —»oo 

(3) 

RUB = S. 


If this is true, we expect that the limit sets R and B will be finite unions 
of angular sectors, as the simulation results shown in Figured suggest. The 
sizes and directions of these angular sectors (and even their number) will, we 
expect, be random, with distributions depending on the initial configuration. 
This is illustrated by simulation results summarized in Figured with initial 
configuration 


R(0) = {(-2,-2), (2,2)}, 
R(0) = {(-2,2), (2,-2)}. 
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Figure 1. Competition Model with -R(O) = {(0,0)}, B( 0) = 
{( 1 , 0 )}. 



Three time progressive snap shots of the process were taken. The plots in 
the figure [3] suggest that stabilization of the shape was taking place on the 
considered time interval. 


2. The Competition Model on Z 1 

The coexistence problem for the Competition Model in one dimension is 
considerably simpler than in higher dimensions. Since the limit shape of the 
Richardson model in one dimension is an interval, no auxiliary hypothesis 
is needed. 

Proposition 1 . For any nondegenerate finite initial configuration on Z, the 
event of mutual survival in the two-species Competition Model has positive 
probability. 


Proof. Without loss of generality, we may assume that the initial configu¬ 
ration consists of a finite interval of red sites with rightmost point —1 and 
a finite interval of blue sites with leftmost point 0, since [a translate of] 
such a configuration may be reached in finite time, with positive probabil¬ 
ity, from any nondegenerate initial configuration. Let Xf and 1} be the left- 
and right-most occupied sites (of either color) at time t, and let Zt be the 
leftmost blue site. Note that as long as X± < Zt < Yt, there will be both red 
and blue sites: all sites to the left of Zt are red, and all sites to the right are 
blue. Each of the processes X t and Y t is a pure jump process, with jumps 
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Figure 2. Competition Model with R( 0) 

{(-2,-2), (2,2)} and £>(0) = {(—2,2), (2,-2)}. 



of size 1 occurring at rate 1; hence, with probability one, as t —► oo, 

Xf/t —* —1 and 
Y t /t — 1. 

The process Z t behaves, up to the time of hrst exit from ( X t ,Y t ), as a 
continuous-time simple nearest-neighbor random walk on the integers. Con¬ 
sequently, there is positive probability that Z t never exits the interval (X t , Y t ). 
But on this event, both species survive. □ 

This simple argument clearly shows what the difficulty in higher dimen¬ 
sions will be: In one dimension, the interface between (connected) red and 
blue clusters is just a point; but in higher dimensions, it will in general be 
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a hypersurface, whose time evolution will necessarily be somewhat compli¬ 
cated. 


3. Preliminaries 

3.1. Graphical Constructions. The Richardson model, the voter model, 
and the two-species competition model all admit graphical contructions us¬ 
ing percolation structures. Such constructions make certain comparison ar¬ 
guments and duality relations transparent. We briefly review the construc¬ 
tion here, primarily to emphasize that the same percolation structure can be 
used to simultaneously build versions of all three processes with all possible 
initial configurations. See, for instance, |3j for further details in the case of 
the Richardson model and the voter model. 

The percolation structure II is an assignment of independent, rate-1 Pois¬ 
son processes to the directed edges xy of the lattice Z d . (For each pair 
{x,y} of neighboring vertices, there are two directed edges xy and yx.) 
Above each vertex x is drawn a timeline, on which are placed marks at the 
occurrence times T? y of the Poisson processes attached to directed edges 
emanating from x ; at each such mark, an arrow is drawn from x to y. A 
directed path through the percolation structure II may travel upward, at 
speed 1, along any timeline, and may (but does not have to) jump across 
any outward-pointing arrow that it encounters. A reverse path is a directed 
path run backward in time: thus, it moves downward along timelines and 
jumps across inward-pointing arrows. A voter-admissible path is a directed 
path that does not pass any inward-pointing arrows. Observe that for each 
vertex y and each time t > 0 there is a unique voter-admissible path be¬ 
ginning at time 0 and terminating at ( y,t ): its reverse path is gotten by 
traveling downward along timelines, starting at ( y,t ), jumping across all 
inward-pointing arrows encountered along the way. 

Richardson Model: A version Z(t ) of the Richardson model with initial 
configuration Z( 0) = £ is obtained by setting Z(t) to be the set of all 
vertices y such that there is a directed path in the percolation structure II 
that begins at (x,0) for some x G £ and ends at (y,t). 

Voter Model: A version ( R(t),B(t )) of the voter model with initial config¬ 
uration R( 0) = (3, B( 0) = (3 C is gotten by defining B(t) = R(t) c and R(t) 
to be the set of all vertices y such that the unique voter-admissible path 
terminating at (y, t) begins at (x, 0) for some x € (3. 

Two-Species Competition Model: Fix an initial configuration R( 0) = £, 
B( 0) = £. Erase all arrows that lie only on paths that begin at points 
(x,0) such that x 0 £ U £; denote the resulting sub-percolation structure 
% f . Define R(t) (respectively, B(t)) to be the set of all vertices y such that 
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there is a voter-admissible path relative to II^ that ends at (y, t ) and starts 
at (x,0) with x G £ (respectively, x G (). 

The graphical construction yields as by-products comparison principles 
for the Richardson, voter, and competition models. First, the set R(t) U 
B(t) of vertices occupied by either Red or Blue particles at time t in the 
competition model coincides with the set Z(t) of occupied vertices in the 
Richardson model when Z{ 0) = R(0)UR(0). Second, if R^(t), B^\t) is the 
voter model with initial configuration i?W(0) = and R(t),B(t ) is the 
competition model with initial configuration R( 0) = £ D , B(0) = ( C £ c , 
then for all t > 0, 

R(t) D R {1) (t). 


3.2. Voter Model: Invasion Times. How long does it take for a Red 
vertex to be overrun by Blue? Clearly, for either the voter model or the 
competition model the answer will depend, at least in part, on how far away 
the nearest Blue vertices are. The comparison principle implies that, for any 
given value p of the distance to the nearest blue vertex, the worst case (for 
either model) is the voter model with initial configuration R(0) = D(0,p) 
and -B(O) = D(0,p) c , where D(x,p ) denotes the disk of radius p centered at 
x (more precisely, its intersection with the lattice Z rf ). 

Lemma 1. Fix /3 G ( 5 , 1), and denote by ( R(t ), B(t)) the state of the voter 
model at time t. There exist constants ci, C 2 >0 (depending on /3) such that 
for all p > 1 and all t G [0, p\, if R{ 0) contains the disk of radius pr centered 
at the vertex x, then 

P{x £ R(t)} < ci exp{—c 2/ o /3 ~ 1/2 }. 

Remark. This holds for any norm on M d , not just the Euclidean norm: in 
particular, it holds for the Richardson norm defined below. The constant 
C 2 > 0 may, of course, depend on the norm. 


Proof. The dual process of the voter model is the coalescing random walk 
(see [3| or 0). Thus, the probability that the vertex x is blue at time t 
coincides with the probability that a continuous-time simple random walker 
started at x at time 0 will land in the set -B(O) at time t. (This is not difficult 
to deduce directly from the graphical construction above: the event 1 G 5(1) 
occurs if and only if the reverse voter-admissible path started at (x, t) will 
terminate at (y, 0) for some y G B( 0); but the reverse voter-admissible path 
is a simple random walk.) Hence, if R( 0) D D(x, p 13 ) then this probability is 
dominated by the probability that the continuous-time simple random walk 
exits the ball D(x,p@) by time p. □ 
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3.3. Richardson Model: Shape Theorem. The first-order asymptotic 
behavior of the Richardson model on the integer lattice T, d is described by 
the Shape Theorem [2j. Denote by Z(t) the set of vertices of T, d that are 
occupied at time t, and by the probability measure describing the law of 
the process given the initial condition Z( 0) = £. For any subset Z C M. d , 
define 

Z = {x £ : dist(x, Z) < 1/2}, 

where dist denotes distance in the L°°— norm on M. d . For any subset Z C M. d 
and any scalar s > 0, let Z/s = {y/s : y £ Z}. 

Theorem 2 (The Shape Theorem). There exists a nonrandom compact 
convex set S C invariant under permutation of and reflection in the 
coordinate hyperplanes, and with non-empty interior, such that for any finite 
initial configuration Z( 0) = ( and any e > 0, with P^ —probability one, 
eventually (i.e., for all sufficiently large s) 

(1 -e) -S C ^2 C (1 + e) -S. 
s 

The exact shape of the limiting set S remains unknown. A simple argu¬ 
ment shows that S is convex, but nobody has succeeded in proving that it 
is strictly convex. Let | • | be the norm on associated with the shape set 
S, that is, for x £ |®| = infjf : x £ S ■ t}. That this is in fact a norm 

follows from the convexity of S. The Shape Theorem is equivalent to the 
statement that the set of occupied sites grows at speed one, relative to the 
norm | • |, in every direction. 

The Richardson model admits a description as a first passage percolation 
model, as follows. To each edge of the lattice Tfl, attach a mean one expo¬ 
nential random variable, the “ passage time”, in such a way that the passage 
times of distinct edges are mutually independent. For any self-avoiding path 
7 , define the traversal time t( 7 ) to be the sum of the passage times of the 
edges in 7 . For any finite set f of vertices and any vertex x, define the pas¬ 
sage time T(f,x) from f to x to be the infimum of the traversal times r( 7 ) 
of all self-avoiding paths connecting 1 to (. A version of the Richardson 
model Z(t) with initial configuration Z(0) = ( is given by 

Z(t) := {x £ Z d : T{(,x) < t}. 

The first-passage percolation representation gives simultaneous realizations 
of Richardson evolutions for all initial configurations. Since traversal times 
of paths are the same backwards and forwards, the following duality property 
is immediate: For any finite subsets F,G C Z d and any t > 0, 


(4) 


p F {z(t)nG = 0} = p G {z{t) n f =0}. 
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Kesten 15] and Alexander JT] have established large deviation results for 
the passage times in first passage percolation that specialize to the Richard¬ 
son model as follows. 

Theorem 3. There exist constants c\ and C 2 > 0 such that for any e > 0, 
y E Z d and \y\ 1/2+e < t < \y\ 3/2 ~ e , 

P{\T(0,y) - |y|| < t} > 1 - Cl • exp{—c 2 f/ \/\y\}- 

3.4. A Ti'iangle Inequality. The hypothesis of Theorem |T| is that the 
Richardson shape S is uniformly curved , that is, that there exists g > 0 
such that for each x E dS there is a (Euclidean) ball D*(x) of radius g 
containing S with x on its surface. Denote by ir : M d \ {0} —► dS the 
natural projection onto the boundary of the Richardson shape, that is, for 
any x E such that i^O, 


7 TX = x/\x\. 

Lemma 2. Suppose that S is uniformly curved, then there exists a constant 
c > 0 such that, for all x S and y E S 

\x — y\ > \x — irx\ + c • | 7 rx — y | 2 


Proof. Let | ■ 11 2 be the Euclidean norm (L 2 -norm) on M f b Any two norms on 
W l are equivalent, and so the Euclidean norm is equivalent to the Richardson 
norm: in particular, there is a constant S > 1 such that, for any x E M d , 

-r • ||x|| 2 < \x\ < 5 ■ ||x|| 2 . 

0 

Let l be the tangent hyperplane to D*(x) at irx. For any y E S, denote by 
y p the (orthogonal) projection of y on l. Since S is uniformly curved (and 
hence strictly convex), 

\x-y\>\x- 7rx\ + |y - y p \. 

Elementary trigonometric observations imply that, since y E D*{x) and y p 
is on the tangent line l, there exists a constant Co that depends only on g 
and such that for all such y’s 

II y - Vph > Co • \\-nx - y\\l 

Hence, 

\y-y P \ >( p- Kx - y\ 2 

which proves the inequality of the lemma. □ 
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4. Proof of Theorem □ 

4.1. Strategy. We begin by showing that it suffices to restrict attention 
to a special class of initial configurations, which we dub sliced Richardson 
shapes. These are obtained as follows: Run the Richardson model (starting 
from the initial configuration with two adjacent occupied sites, one at the 
origin) for a (large) time T, and let Zt be the occupied set. Let be the 
subset of Zt consisting of all points with positive first coordinates, and let 
Qt = Zt \ £. Observe that, starting from any nondegenerate finite initial 
configuration the competition model can evolve to a sliced Richardson shape 
in finite time, with positive probability. (This will occur if, following the first 
time that there are adjacent Red and Blue sites, only these sites reproduce, 
and only on their sides of the hyperplane separating them.) Thus, it suffices 
to prove that for all sufficiently large T, with positive probability the sliced 
Richardson shape Ct is such that 

(5) P^{n t >o-R(t) = 0 or n t > 0 B(t) = 0 } < 4e. 

Here and in the sequel P^^ will denote the probability measure governing 
the evolution of the competition model under the initial condition R( 0) = £, 
-B(O) = C- By the Bonferroni inequality, it suffices to prove that 

( 6 ) P^{n t > 0 R(t) = 0 } < 2e. 

The idea behind the proof of © is this: If the initial condition £, ( is such 
that £ and ( are, approximately, the intersections of kS with complementary 
angular sectors A, A c in based at the origin, for large k , then at time 
t = 5k the sets R(t),B(t) should, with high probability, be approximately 
the intersections of k( 1 + 5)S with the same angular sectors A, A c . This 
is because (1) the Shape Theorem for the Richardson model implies that 
R(t) U B(t) should be close to (1 + 5)kS; (2) the uniform curvature of S 
implies that the first occupations of vertices in A D ((f + k)S \ kS ) and 
A c n ((t + k)S \ kS) should (except for those near the boundaries) be by 
Red and Blue, respectively; and (3) Lemma [T] implies that, once a region is 
totally occupied by Red, it must remain so (except near its boundary) for a 
substantial amount of time afterward. 

4.2. Stabilization Estimate. The key step is to show that once one of 
the species (say Red) has occupied an angular sector in the Richardson 
shape, it is very unlikely for the opposite species (Blue) to make a large 
incursion into this sector for some time afterward. Henceforth, let d be 
the metric associated with the Richardson norm. For any set ( and any 
vertex x, define the distance d(£, y ) between ( and x to be the infimum of 
the distances d(y , x) for all vertices in y £ £. For any point x £ and 
any r > 0, denote by D(x\r) the disk of radius r centered at x relative 
to the metric d. (We shall not attempt to distinguish between open and 
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closed disks, as this distinction will not matter in any of the estimates.) For 
?T < i "2 denote by D(x;r\,r 2 ) the annular region D(x;r 2) \ D(x;ri). For 
each 2 G dS and any g > 0, dehne the angular sector A(z; g) of aperture g 
centered at z by 

A(z; g) := {y G R d \ {0} : d(ny, z) < £>}. 

Fix r > 0, n > 1, /3 G (1/2,1) and a G (1/2,1) such that (/ 3 + l)/2 < a , 
and let A\ C A 2 be angular sectors with common center z and apertures 
r < r + n a_1 , respectively. Fix 8 G (0,1), and define 

7 l 0 = 1Z% = D( 0; n/( 1 + <5), n - n^) fld 2 n Z d , 

B 0 = Bf = (D( 0; n/(l + 5)) U {D( 0; n/(l + 8),n + vP) D A c 2 )) n Z d , 

Bi = Bf = (D(0;n)Ud5)nZ d , 

n x = TCl = {D{ 0; n, n(l + <5) - (n + dn)^) D Ai) D Z d , 

Lemma 3. There exist constants ci,C 2 > 0 such that the following is true, 
for any n > 1. If the initial configuration £, £ is such that f D 77q and 
( C BJ, then 

(7) 1 — P^^{B(8n) C Bf} < cin M exp{—C 2 (dn) /3-1 / 2 }. 

Proof. To prove 0 we find exponential upper bounds on 
Pn 0 ,B 0 i B {Sn) n 77i / 0} and 
P^BoiBiSn) n(77iU Si)'V 01- 

Claim 1. For all sufficiently large n and for all x G 77 1 

d(x, 77 q) + 4(<5n)^ < d(x, Bo) 


Proof. First, observe that by Lemma El for every x G 1Z\ such that \x\ > 
n + , we have 

d(x, Bo) > \x — nx ■ (n + nP)\ + c ■ n 2a_1 . 


Also, 


d(x, 77o) < | a: — ttx • (n + nP )| + 2 jp. 

These inequalities imply Claim |T| for all x G 77 1 D D(0;n + n^) c and all 
sufficiently large n. 


Next, if x G 77i is such that n < \x\ < n + n^, then 

d(x, IZo) < Sn/ 3 . 

Also, by Lemma El for all y G Bo, 

d(y , x) > c ■ (8n] 2a ~ l . 

These two inequalities imply Claim El for all x G 77 1 n 77(0; n,n + n@) and 
all sufficiently large n. □ 
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Claim 2. With probability —► 1 as n —> oo, /or every x E T^i, every site 
of the ball D(x; (Sri)*) will be colonized by time t x = d(x,H o) + 2 (Sri)*. In 
particular, there exist constants c\ > 0 and c 2 > 0 (not depending on x) 
such that for every x E H\ 

Pn 0 ,B 0 { D (x, ( Sn ) /3 ) £ ^(tz) U .B(t x )} < ci ■ (<5n ) d/3 ■ exp{-c 2 • (Sn)^ 1/2 }. 
Proof. Fix x E H\ and let t x = d(x, Ho)+2(Sn) 33 . For every z E D(x; (Sn) 13 ), 
(8) d(z, Ho) + (Sn) 13 < t x . 

Notice that the number of sites in D(x\ (Sn) 33 ) is of order (5n) d33 . By Theo¬ 
rem ISJand ©, it follows that for some ci > 0 and c 2 > 0 (not depending on 
x) 

P-r 0 ,b 0 { d ( x 'i (Sn) 13 ) <£_ R(t x ) U B(t x )} = Pn 0 uB 0 { D ( x ; (Sn) 33 ) (jf Z(t x )} < 

< ci • (Sn) d33 exp{—c 2 • (Sn) 33 ^ 132 }, 

Hence, with probability —» 1 as n —» oo, for every x E H\, every site of the 
ball D(x\ (Sri) 33 ) will be colonized at time t x . □ 

Define the boundary dA of a set A C Z d as the set of all z E A that have 
at least one nearest neighbor that is not in A. Next, for a set A C Z d let 
t(A) be the first time at which the blue species reaches A. 

Claim 3. With probability —> 1 as n —> oo, for every x E H\, the blue 
species will not reach the ball D(x](Sn) 33 ) by time t x . In particular, there 
exist constants c\ > 0 and c 2 > 0 such that for every x E Hi 

P-r,o,B 0 {t(D(x-, (Sn) 33 )) < t x } < ci ■ (Jn) (d_1)(/3+1) • exp{-c 2 • (i Sn ) /3_1 / 2 }. 
Proof. Notice that 

Pn 0 ,B 0 {r(D(x- (Sn) 0 )) < t x } < P Bo {Z(t x ) n D(x\ (Sn) 33 ) / 0} 

= Pd(x;(5ti)P){Z( T x) H Bo A 0} 

By Claim El for large n we have d(x,Bo) > t x + 2 (Sn)*. Hence, 

P D(x-,(Sn)f>){ Z i T x) n B 0 + 0} < P D (x\{5n)P){Z( T x) <t D(x;T x + 2(Sn) 33 )}. 
Obviously, for every z E dD(x ; (Sn) 33 ) and y E dD(x;r x + 2(Sn) 33 ), we have 

d(z,y) >t x + (Sn) 33 /2. 

Apply Theorem El to each pair of such vertices to get: 

P{T(z,y) < t x } < ci ■ exp{—c 2 • ( Sn) 33 ~ 1/2 }. 

The number of vertices in dD(x, (Sn) 33 ) is of order (Sn) 3d ~ 1333 , and the num¬ 
ber of vertices in dD(x,r x + 2 (Sri) 33 ) is of order at most (Sn) 3d ~~ 13 . Hence, 

P D(x-,(Sn)^{ z (n) t D(x-,r x +2(5n) 33 ) < cv(Sn) 3d - l) ^ +1) exp{-c 2 -(J?r)^ 1/2 }. 
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This finishes the proof of Claim 01 □ 

Claim 4. With probability —> 1 as n —> oo, there are no Blue particles in 
1Zi at time 6n. In particular, there exist constants c\ > 0 and C 2 > 0 such 
that 

Piz 0 , Sol'll n B(Sn ) / 0} < ci • (dn)( d ~ 1 ^ /3+1 ) +d • exp{-C 2 • (dn)^ 1 / 2 }. 

Proof. By Claim 0 for all x € TZ\ with t x > Sn, 

(9) Ptz 0 ,b 0 {x e B{Sn)} < a ■ (<5n) (o! ~ 1)(/3+1) • exp{-c 2 ■ {Sn)P~ l/2 }. 

Next, for all x € IZi with t x < Sn, Claim[2]and Claim [3] imply that for some 
ci > 0 and C 2 > 0 

PtZo,b 0 {D(x, ( Sn ) /3 ) <£ R(t x )} < ci • (<5n) (d_1)(/3+1) ■ exp{-c 2 • (■ Sn ) /3_1/2 }. 

Hence, by Lemma El there exist constants ci and C 2 such that for every such 
x we have 

(10) Pji 0 ,B 0 {x 0 R(Sn)} < ci • (<5n) (d_1)(/3+1) • exp{-c 2 • (< Sn)^~ 1/2 }. 

The number of vertices in 1Z\ is of order ( Sn) d . Thus, by combining Q and 
m, we get 

p Ko,Bo{ B ( Sn ) n III / 0} < Cl • (<5n) (d_1)(/3+1)+<i ■ exp{—c 2 ■ (Sn) l3 ^ 1/2 }. 

□ 

Claim 5. With probability —* 1 as n —* oo, the blue species will not reach 
the set (TZ\ U £>i) c by time Sn. In particular, there exist constants c\ > 0 
and C 2 > 0 such that 

Pn 0) B 0 {x{(Ili u ^i) C ) ^ <W ^ c i ' ™ 2(d_1) • exp{-c 2 • (Sn)P~ 1/2 }. 

Proof. For large n the distance between the sets Bq and (IZ\ L)Bi) c is greater 
than Sn + (Sn)@. The number of vertices on the boundary of Bo is of order 
n^ 1 . Using the same line of argument as in the proof of ClaireElwe get 

PnoPoMCRi U £>i) c ) < Sn} < P Bo {Z(Sn ) <£K i U Hi} < 

ci • n 2 ^ -1 ) • exp{—C 2 ■ (Sn)P~ 1 / 2 }. 

□ 


Now, Claim El and Claim El imply 0 and finish the proof of Lemma El 
PTlo,Bo{ B (5n) <£ Pi} < ci ■ n 3d • exp{—c 2 • ( Sn ) /3_1/2 }. 


□ 
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4.3. Proof of ®. Let Z t be the set of sites occupied by the Richardson 
evolution (started from the default initial configuration) at time t. Fix 
T = Tq > 1 and 5 > 0, and set 


T n = T{ 1 + <5)-, 

t n = T5(l + 5) n ~ 1 , and 


= r n — t = y] t r 

3 = 1 


Fix f3 £ (1/2,1), and for each n = 0,1,... define events 

K : = {(1 - Tt 1 ) • 5 C Z(T n )/T n C (1 + Tf _1 ) • 5}, 
G n : = {(!- Tt 1 ) • 5 C Z(r n )/T n C (1 + T^ 1 ) • 5}. 


By the Kesten-Alexander large deviation theorems (Theorem 0 , 

OO 

T i™ E p «) = ° 

n=0 

and so, for sufficiently large T, the probability is nearly 1 that the configu¬ 
ration Zt will be such that 

OO 

( 11 ) Y, P ^Gn)<e. 

n=0 


Fix an initial configuration Z so that the preceding estimate holds for 
Zt = Z, and use this to construct the split Richardson shape £ = £r, C = Ct 
as in section rm above: £ and £ are the subsets of Zt with positive and 
nonpositive first coordinates, respectively. Since the union of the Red and 
Blue sites in the competition model evolves as the Richardson model, it 
follows from EH) that if the initial configuration is R{ 0) = £, B(0) = £ then 
with probability in excess of 1 — e, 

(12) (1 - T%~ 1 ) ■ S c C (1 + rf- 1 ) • 5 

n 

for all n = 0,1, 2,... . Denote by G* the event that m holds for all n > 0. 

On the event G*, the union of the Red and Blue regions will, at each time 
T n — T, fill a region close enough to a Richardson shape that the estimate 0 
will be applicable whenever the Red and Blue populations are restricted (at 
least approximately) to angular sectors. Thus, define sequences A™ C A% of 
concentric angular sectors with apertures r\{n) < r 2 (n) such that 

r 2 (n) — ri(n) = and 

r 2 (n + 1) = n(n) 

and with r 2 ( 0 ) chosen so that A® is the halfspace consisting of all points 
in with positive first coordinates. Here ^</3<(/3 + l)/2<a<l 
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as in Section ro above. Note that the second equality guarantees that 
A% +1 = A™. This in turn, together with the fact that the sequence T n is 
increasing, implies that that the angular sectors are nested: A™ +1 C A™. 
Moreover, because T n is an exponentially growing sequence and a < 1, 

lim rr(n) := > 0 

n —xx) 

provided T = To is sufficiently large. Therefore, the intersection 

,4°° = n ™ =1 A™ 

is an angular sector with nonempty interior. 


Finally, for each n > 0 define H n to be the event that at time r n there 
are no blue sites in A% outside the (Richardson norm) disk of radius T n _ i. 
(For n = 0, set T_i = To/(l + <5).) On the event G n (~l H n , the set of all 
occupied sites is close to T n ■ S, and the red sites fill at least the outer layer 
of this set in the sector A%. We claim that for all sufficiently large T, 


(13) 


B 


Wr f)(G n nH n ) 


\n=0 


> 1 - 2e. 


To see this, let v be the smallest index n such that (G(( U H (() occurs. Since 
v = n can only occur on G n -\ n H n _ i, 


{ ,y — n l < I G n -i PI H n _ i). 

Inequality m provides a bound on the sum of the first of these terms, and 
Lemma |5] bounds the second. Thus, for T sufficiently large, 

OO 

n =0 

this proves m. 


On the event G n (~l H ni the Red species must at time r n occupy at least 
the outer layer of the occupied set in the angular sector A°°. Consequently, 
on the event H n >i (G n H H n ), Red survives! This proves ©■ □ 


5. Concluding Remark 

The preceding argument, in addition to proving that the event that mu¬ 
tual survival has positive probability, also goes part of the way towards 
proving Conjecture Q If at a large time T one of the colors (say Red) occu¬ 
pies the outer layer of an angular sector, then with conditional probability 
approaching 1 as T —> oo it will occupy a slightly smaller angular sector 
forever after. Since the same is true for the other species, it follows that in 
at least some evolutions Red and Blue will each occupy angular sectors. 
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Unfortunately, it remains unclear what happens near the interface at large 
times. Although the preceding arguments show that neither Red nor Blue 
can make too deep an incursion into the other species’ sector(s), it may 
be possible for one to repeatedly make small incursions across the interface 
that engender more (and necessarily thinner) angular sectors in its zone 
of occupation. Thus, it may be that the limit shapes exist, but consist of 
countably many angular sectors. 

Finally, it remains unclear if stabilization must eventually occur on the 
event of mutual survival, that is, if it is necessarily the case that at large 
times T the outer layer of the occupied region must segregate into well- 
defined Red and Blue zones. Since local coalescence occurs in the voter 
model, one naturally expects that the same will be true in the competition 
model; thus far, we have been unable to prove this. 
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